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The operator product expansion is used to compute the matrix elements of composite renormalized operators 
on the lattice. We study the product of two fundamental fields in the two-dimensional cr-model and discuss the 
possible sources of systematic errors. The key problem turns out to be the violation of asymptotic scaling. 



The Operator Product Expansion 
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is widely thought to hold beyond perturbation 
theory. 

The use of Eq. (|l]) in lattice simulations [0-^ 
is still in its infancy. In Ref. |l[] it was suggested 
to use Eq. (|l|) in order to compute renormalized 
matrix elements. The OPE approach consists in 
the following steps: one computes the (matrix el- 
ement of the) l.h.s. of Eq. ([^), then rcnormalizes 
A and B in some scheme, and finally obtains (the 
matrix element of) C through a fit, using some 
perturbative approximation of the Wilson coeffi- 
cients. 

The main sources of systematic errors in this 
approach are the following: 

(a) , finite-size effects and corrections to scaling, 

i.e. lattice artifacts; 

(b) . "power-correction effects" which are due to 

the fact that we truncate the expansion (|^) 
to some finite order in x^: 



(c). 



corrections to asymptotic scaling which 
must be taken in account since the Wilson 
coefficients in Eq. (^ have to be substi- 
tuted by the first few terms of their pertur- 
bative expansion. 



need more explanations. Here we shall focus on 
errors of type |(b)| and |(c)| . 

The use of Eq. (|l|) in a continuum scheme, 
for which only a perturbative computation of the 
Wilson coefficients is available, poses a restriction 
on the operators which can be obtained in this ap- 
proach. Only the operators of lowest dimension 
for each spin sector, i.e. those of lowest twist, 
can be computed using Eq. (^. Higher- twist 
operators give rise to systematic errors of order 
0{x^). Moreover, because of statistical errors, 
only the operators of low dimension can be reli- 
ably extracted: the remaining ones are strongly 
subleading in the region of validity of the OPE. 

Problem [c) can be stated in a cleaner way if 
we get rid of the scale dependence which is in- 
troduced in this approach somehow artificially 
through the Wilson coefficients. One can rewrite 
Eq. (|l|) by making use of renormalization-group 
invariant operators defined as follows: 
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The Wilson coefficients obviously become fi- 
independent and their general perturbative form 
is 
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Errors of type (a) are widely studied and do not 



where A is the "lambda parameter" of the theory. 
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The use of a truncation of Eq. M) introduces 
systematic errors of order 0(log~ (Aa;)). No- 
tice, however, that this approach aUows to com- 
pute directly "infinite-energy" quantities (i.e. the 
renormahzation-group invariant matrix elements) 
which are of interest in phenomenological appli- 
cations. Errors of order log~'^(Aa;) arise also in 
the widely used "non-perturbative renormaliza- 
tion method" 0| in which perturbation theory 
is used to "evolve" the renormalization constants 
computed at some energy scale achievable on the 
lattice up to high energies. 

We have considered several products of opera- 
tors for the 0{N) nonlinear a- model in two di- 
mensions, with lattice action 

S{<r)^^J2(9,cT)l, (5) 

a;, /J 

where a e 5'^"^ = fx+^-fx, and = 3. 

Here we shall refer to the following (respectively 
scalar and symmetric) products of fundamental 
fields: 



(t{:x)-(t{-x)^Wq{x)+0{x^) , 
(7'^{x)(7^{-x) + [a ^b) ~ trace ^ 
r^W2{x) [cr'^cr'' - trace] (0) -I 



0{x^ 
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The Monte Carlo data presented refer to two lat- 
tices: the first one of size L x T ^ 128 x 256 and 
correlation length {am)~^ = 13.632(6) (m is the 
mass gap); the second with L x T = 256 x 512 
and correlation length (am)~^ — 27.094(43). 

The expectation value of the product (|^) be- 
tween states of momentum p is shown in Fig. |l| 
for the two different lattices: corrections to scal- 
ing, i.e. errors of type (a), are completely under 



control in our simulations. 

Our general procedure consists in choosing a 
truncation of the expansion (|^) and in using it 
to fit the data in the region p < \x\ < R. The 
results are independent of p for 1.5 a < p < 3 a. 
We use the stability of the fit with respect to the 
truncation and to the choice of i? as a criterion 
to distinguish whether errors of type (_b] and (c] 
are relevant or not. 

The quality of the fits obtained is well repre- 
sented by the curves shown in Fig. ^ where the 
two-loop expression was used for the first term in 
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Figure 1. The product (0) on the slice xq = 0. 
Circles are the data at (ma)~^ ~ 13.632 and stars 
at (ma)~^ ~ 27.094 (properly rescaled). Contin- 
uous and dashed lines are the fitting curves re- 
spectively with and without 0{x'^) terms. 

the expansion (|^) and the tree-level form for the 
terms of order 0{x^). 

An additional conclusion can be drawn from 
in order to describe the symmetric prod- 

it is ncc- 



Fig. 

uct (|7) up to distances 2x ^ m 
essary (and almost sufficient) to include terms of 
order 0{x^). However this does not mean that 
the terms of order x'^ with the same symmetry of 
the leading one — the higher-twist terms — can 
be obtained from the fit. Indeed, their matrix el- 
ements extracted from the fits are very unstable 
with respect to changes of R. 

In order to test the stability of the fit, we 
considered the vacuum expectation value of the 
product (|^), that is the two-point function. Here 
the renormalization constant of the field is a fit 
parameter. The results are reported in Fig. ^ 
and refer to the lattice with correlation length 
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Figure 2. The field-renormalization constant as 
computed with different truncations of the OPE. 
From top to bottom: one loop, two loops, three 
loops, three loops plus 0{x^) corrections. 

{ma)~^ ~ 27.094; we used /la — 10. As a manifes- 
tation of asymptotic freedom, the various curves 
shrink when R 0. Their i?-dependence be- 
comes weaker as more terms of the perturba- 
tive expansion are included. Nevertheless, even 
if we use the three-loop Wilson coefficient, we 
do not obtain a value of Z independent of R 
within the statistical errors, as it should be in 
the asymptotic-scaling regime. Two remarks are 
in order here: first, it is notoriously difficult to 
reach asymptotic scaling in the 0(3) nonlinear cr- 
model second, we are studing a small effect 
which could be negligible in QCD applications 
with respect to other errors. 

Note that the addition of 0{x^) terms, which, 
in this case, have spin and are therefore higher 
twists, seems to improve the situation, making 
the curve flatter. However, this must be re- 
garded as a spurious effect. Higher twists are 
simply mimicking the contribution of higher or- 
ders in perturbation theory: otherwise, the curves 
obtained without including them should show a 
(mi?)^ behavior. From Fig. |2| we can estimate 
the systematic error due to the use of two- or 
three-loop Wilson coefficients to fit the data in 
this range of mR: the error is approximately 5%. 
We report in Fig. ^ the results of the fit pre- 
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Figure 3. The expectation value {p\a°'a^—tTa,ce\p) 
with momentum: p = (circles); p = 2'k/L (filled 
diamonds); p = Air/L (triangles). Two different 
fits are reported: 0{x'^) corrections are not con- 
sidered in fit (A), while they are included in fit 
(B). 

sented in Fig. ^ Each type of truncation, in- 
cluding or not 0{x^) corrections, gives a reason- 
able, that is p independent and (almost) R in- 
dependent, answer. The difference between them 
should be interpreted as a violation of asymptotic 
scaling and, indeed, it is of the same magnitude 
of the systematic error estimated for the renor- 
malization constant of the field. 

We thank R. Petronzio for useful discussions. 
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